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INPUT ADMITTANCE OF A COAXIAL TRANSMISSION LINE OPENING ONTO 
A FLAT, DIELECTRIC- COVERED GROUND PLANE 

By Calvin T. Swift 
Langley Research Center 

SUMMARY 

The input admittance of a coaxial transmission line opening onto a flat, dielectric- 
covered ground plane is derived, and computational results are given as a function of 
slab thickness in wavelengths. A dielectric constant of 2.57 was selected for the slab, 
which physically corresponds to that for lucite or polystyrene. 

Surface wave contributions are considered in the calculations, and the results show 
that over some ranges of slab thickness 90 percent of the power supplied to the antenna 
is trapped in the dielectric. 

Comments pertinent to a physically realizable antenna are included. 

INTRODUCTION 

In order to compute the input impedance of an antenna, a difficult electromagnetic 
boundary value problem must be solved. From a practical viewpoint, the problem should 
be solved either in rectangular, cylindrical, or spherical coordinates. This therefore 
requires a simple aperture and feed which can easily be described by these coordinates. 
The simplest example is that considered herein; that is, the large coaxial transmission 
line opening onto a flat ground plane, which can readily be described in cylindrical 
coordinates. 

This antenna was originally investigated by Levine and Papas (ref. 1) who consid- 
ered the same problem, with no dielectric cover. A variational expression for the 
admittance was derived which included the effects of higher order transmission-line 
modes on the admittance. Only the dominant mode was considered in the calculations, 
and the results compared favorably with the measured values of Hartig (ref. 2). 

Later Galejs (ref. 3) treated the annular slot on a flat ground plane, feeding into a 
large circular waveguide containing a stratified plasma. By choosing a waveguide radius 
large compared to a wavelength, Galejs was able to infer the admittance of the coated 
slot on an infinitely large ground plane. 



In this paper, as in reference 1, a flat ground plane of infinite area is initially 
assumed. The vector potential is written in terms of Hankel transforms and a three- 
region boundary value problem is then solved by assuming the existence of only the 
dominant transverse electromagnetic (TEM) mode at the aperture. Parseval's theorem 
is developed, which leads to an admittance expression containing an integral over the 
real axis in mode space, which must be numerically evaluated. Pole contributions in the 
integral are evaluated which physically describe surface wave propagation in the loss- 
free medium considered herein. 

Numerical results are given as a function of dielectric slab thickness in wave- 
lengths for a slab whose dielectric constant is 2.57. 

SYMBOLS AND ABBREVIATIONS 

Rationalized mks units are used for all calculations. 

A* magnetic vector potential 

a radius of inner conductor of transmission line 

b radius of outer conductor of transmission line or normalized input 

susceptance 

E electric field intensity 

g r normalized radiation conductance 

g s normalized surface wave conductance 

H magnetic field intensity 

k = Nk Q 

kg free- space wave number 

k z axial wave number 

N index of refraction of slab 
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p* 

Rn(p) 

t 

x,y,z 

Y 0 

Y n 

^in 

Zc 

z 0 

0 

r 

e 

e 0 

x o 

M 0 


index of refraction of dielectric in coaxial transmission line 
complex conjugate of power 

general radial solution in the waveguide, J 1 (>-nP)Yo(^ n a^ - Jg(A n a )Yi(;\- n p) 
where n is an integer = 1 

time 

Cartesian coordinates 
characteristic admittance of free space, 

]Y 0 N’ 2 

admittance of nth transmission-line mode, 
normalized input admittance 

"I k 

characteristic impedance of coaxial transmission line, — - In — 

’ 2irY 0 N’ a 

slab thickness 

normalized radial mode number, ^/kg 
location of surface wave pole in mode space 
reflection coefficient 
permittivity 

permittivity of free space 
wavelength in free space 
cut-off wave number for higher order modes 
wavelength in dielectric slab 
permeability of free space 
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radial wave number 


k 

£ radial wave number used as dummy variable 

p,cp , z cylindrical coordinates 

co angular frequency 

Subscripts: 

p,cp,z principal vector components in cylindrical coordinates 

ap aperture 

Superscripts: 

0 within coaxial line 

1 in dielectric slab 

II in free space 

* complex conjugate except for vector potential 

Abbreviations: 

TEM transverse electromagnetic 

TM transverse magnetic 

VSWR voltage standing-wave ratio 

A bar over a symbol indicates the Hankel transform. The letters PP before an 
integral indicate the principal part of the integral. 

THEORY 

A coaxial transmission line of inner radius a and outer radius b, as shown in 
figure 1, opens onto a conducting ground plane coated with dielectric material of 
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thickness Zq. The line is excited in the dominant TEM mode and a time-harmonic field 
of the form e - ^ w * is assumed. Since the fields are symmetric in cp and TEM to z 
inside the guide, the only nonvanishing field components are H^, E p, and E z . As such, 
the fields may be developed from a magnetic vector potential consisting of only a 
(^-component. From Stratton (pp. 23-34 of ref. 4), this potential A^O^z) can be related 
to the field components as follows: 


H<p(p,z) = jwA* (p,z) 

1 9Aj(p,z) 


En(p,z) — 


9z 


(i) 


Ez(p,z) = -“^[pAj(p,z)] 


where, for source-free homogeneous media, A* (p,z) satisfies the Helmholtz equation 




8p‘ 


+ ~ 


P 9p 


8z* 


+ k' 


A* =0 


2 <P 


( 2 ) 


The structure is unbounded in the radial direction; therefore, the solution to equation (2) 
is: 

Aj(p,z) = (3) 

where A^(£,z) is the Hankel transform of A^(p,z), that is, 

A*(£,z) = j^ pA*(p,z)J 1 (|p)dp (4) 

The substitution of equation (3) into equation (2) reduces the Helmholtz equation to 

= 0 < 5 > 

In free space (z g z^j, only outgoing traveling waves exist; therefore, this solution 
to equation (5) is 

—II n 

A* (£,z) = C^e 3 z (6) 

where the roots of k z ^ must be chosen so that 
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k z n = k 0 f - («Ao) 2 

(l«Ao| « ‘ 

) pa) 

k z“ = i^Aof - 1 

(|«/ k o| > \ 

) (1b) 

Standing waves are supported within the homogeneous dielectric slab (0 £ z g z^ 
this solution to equation (5) is 

j); hence, 

aJ^z) = A(£)e ]kz Z + BU)e~ 

i k z Iz 

(8) 

where, for real index of refraction N, the branch of k z * 

is arbitrarily selected so that 

k z‘ = 2 - (V k o) 2 

^Trr 

o 

A 

(9a) 

kz 1 - Jko^/ko) 2 - 1,2 

(l*A)| >N ) 

(9b) 


If Ep (£) is defined as the unknown transform of the electric field at z = 0, 
the coefficients A(£), B(§), and C(£) may be eliminated in preference to E p ap (£) by 
means of the boundary conditions at z = 0 and z = z Q ; that is, 


V t («> z o)"V( { ' Z o) 

E p n (^o) = V(«> z o) > < 10 > 

VpHm = i p>ap ({) 


By using equations (10) to solve for A(|), B(£), and C(|), the transforms of the 
magnetic field in regions I and n may be obtained in the following form: 


V(£, Z ) = 


_ kp N ^ e o/^Q E p,ap^^z^ cos k z ( z " z o) + i N k z sin k z ( z " z q 1 
k z ^ N 2 k z II cos kjj^Q - jk z *sin k z I zJJ 


s n,., (Z " 0) 

’ 2 II I I I 

N k z cos k z z Q - jk z sin k z z Q 
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The general solution of the fields inside the transmission line (z < 0) is a super- 
position of a forward propagating TEM mode plus all higher order of TM reflecting 
modes; that is, 



(12a) 

(12b) 


where B^ip) is the radial solution and X n is the cut-off wave number of the nth 
higher order mode, as determined from the transcendental equation 

J 0 ( X n a )Y 0 (*nb) - Y Q (x n a)J 0 (x n b) = 0 (13) 

which follows from the requirement that E z must vanish at p = a and p = b. 

From the continuity of Ep and By within the aperture a rigorous solution can 
be obtained for the unknowns Ep ap(£), E, and A n . However, for the present computa- 

2 ? 

tions, it will be assumed that the quantity A n - k is large enough so that the sum of 
all higher order mode contributions is negligible (Aj = A 2 = . . . = A n = 0). Under this 
assumption, then, the internal fields at the aperture are 

E p°(p,°) = ^(1 + T) (14a) 

n N ’ Y 0 A 0 

H^°(p,0) = - r) (14b) 


where r is the reflection coefficient. 

Taking the transform of equation (14a) immediately yields 

r b [JoUfc) - Jo(£a)l 

Ep,ap(D = A q (1 + T)j JiUp)dp = -A 0 (l + F) b (15) 

The complex conjugate of power flow out of the transmission line at z = 0, in terms of 
the field quantities in region 0, is given by 

P * = I J ^ Ep°*(p,0)H (p °(p,0)p dp dcp = tt|A 0 | 2 (1 + r*)(l - T)Y 0 N' In | (16) 

0 a 
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And, from Parseval's theorem (appendix A), the complex conjugate of power flow into 
region I, in terms of the field quantities in region I, is given by 




poo |Eq ap(^)| kz cos k z zq - jN k z sin kg zg £ 

= Trk 0 N 2 Y 0 L (17) 

° k z I jN 2 k z II cos k z I z Q - jkjsin k^z^ 

Substituting explicitly for E n „ n (£) (eq. (15)), and requiring continuity of power across 
the aperture results in 


'in 


1 - r 
i +r 


n 2 r " 

N' In S J 0 
a 


1 . , N 2 Jl.- ^ 2 tan 

\|n 2 - |3 2 

(kgZO^ 2 - f 2 } 

Jjg (kgb/3j - J Q ^k Q a/3| d/3 

jS^N 2 - /S 2 

VN 2 - /3 2 


(18) 


where /3 = |pQ. Equation (18) is valid for complex N. For N real and greater than 
one, it is convenient to break up the integral in equation (18) as follows: 


yin = 


N 2 


N* In 

b 

a 


p 1 pN p 

+ F 2 (i3)di3 + J N F 3 (/3)d)3 


(19) 


The choice of the phases of k * and k ^ (eqs. (9) and (7), respectively) requires that 

z z 

1 - i N ^I^ tan(vo'f>^|[jo(ty>0) - J o( k O a ^)] 2 

\j hr - 0* 


F 2 (/3) = 


P'Jn 


2_ fi 2 


^^-itan( Vo ^T7j 


_/n 2 - /3 2 


(20a) 
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For the case of no coating ^Zq = 0) 

n 2 f ((3) . [ J oM ' J oM)l 

1 (3'/W 2 


(20b) 


(20c) 


(21a) 


and 

N 2 F 2 (/3) = N 2 F 3 (i3) = 

which checks with the results of Levine and Papas (ref. 1). 

The contour of integration of the integral (eq. (19)) is shown in figure 2. In addi- 
tion to the branch point at /3 = 1, a finite number of poles exist on the real axis between 
j3 = 1 and /3 = N as a result of nonconvergence of the Hankel transform for A^,(p,z). 
These poles are of order one, and they occur at the values of /3 for which the denomina- 
tor of equation (20b) vanishes. These poles physically describe the launching of surface 
waves (i.e., energy which is confined within the dielectric coating). This represents 
real power and may be computed from tt j times the residue of equation (20b). Thus 
the complete expression for the input admittance is 

y in = N2 ~ f 1 F,(/3)d/3 + PP f N F 2 (/3)d|3 + f F 3 (/8)d/3 + irj Y Res = g r - jb + g s (22) 

N' ln^ J 0 J 1 J N 6 

a u J 


-{ j oM- j oM] 


m* - 1 


where g r is the normalized radiation conductance, g s is the normalized surface wave 
conductance, and b is the normalized susceptance. The PP preceding the integral of 
F 2 (0) denotes the principal part. 
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From appendix B, the expression for the surface wave conductance is given by 


S s 


W 


N'ln| 


7rj / Res = 


7tN 2 V 

(Vo) 1 *’ ln |n 




[ J o( k Q b ^n) - J 0 ( k 0 a ^n)] 2 


1 + 


'n 2 - 1 ) 

sin 

2k 0 z 0 

(n 2 - fi„ 2 ) 

- V 

( 2k o z o\f 

N 2 - fin 2 ) 



where /3 n is the nth root of the denominator of equation (20b); that is, 



(23) 


(24) 


• It is to be noted that an onset of a surface wave mode occurs where /3 n = 1 and 
the argument of the tangent is 0, v, 2v } . . that is, where the thickness in wave- 
lengths satisfies the equation 


z 0 _ n 

X 0 2\/n S - 1 


(n = 0, 1, 2, . . .) (25) 


which is the same cut-off condition for TM plane waves propagating in a grounded dielec- 
tric sheet. It is of interest to note that the surface wave conductance becomes smaller 
as kpb and k^a become larger. This would be expected, because the radiation 
becomes more directive as the aperture size becomes large compared to a wavelength. 

If the dielectric slab contains a small amount of loss, the poles will be displaced slightly 
off the real axis, and thereby will give no residue contribution. The complex integrand 
of equation (18) will, however, be large in the vicinity of the poles and when integrated 
should give a conductance contribution close to that of g s for the lossless slab. 

As the losses increase, the poles and branch point will be far removed from the 
real axis, and the integrand of equation (18) will be a smooth function, presenting no 
problem for numerical evaluation. 


RESULTS 

Equation (22) was numerically evaluated on an IBM 7094 electronic data processing 
system. The Gaussian quadrature method was used to evaluate the integrals, and the 
surface wave conductance was evaluated separately by using equations (23) and (24). 

A dielectric constant of 2.57 was selected for the exterior slab and 2.00 was 
selected for the material inside the line. The ratio of the radius of the outer conductor 
to the inner one b/a was taken to be 2.00. Nine values of kga were chosen in the 
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range 0.60 ^ k^a S 2.00. For each value of k^a, admittance computations were per- 
formed as a function of z q/X 6 , where Zq is the thickness of the slab and X € is the 

z n i n 

wavelength in the dielectric slab. In the region 0?^ = the incremental values of 

Ag OZ 

3 Zo 1 7 

1/32 and 1/16 in the interval — = Two additional calculations were 

lb 


were 


Z 0/ X e 

zn ^ po 

performed at ^ = — and One surface wave pole exists for 0 ^ ^ ^ 0.64 and two 


z 0 
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poles exist for the remaining computations. 


Plots of the input admittance as a function of slab thickness in wavelengths are 
shown in figure 3. For the smaller values of kg a, the radiation conductance is small 
and relatively constant over the entire range of dielectric thickness considered. How- 
ever, as kpa increases, the radiation conductance increases considerably and begins 
to vary periodically with slab thickness. It is of interest to note that, for thicknesses 
beyond a quarter-wavelength, maximums and minimums occur at roughly every quarter- 
wavelength of slab thickness. 


The surface wave conductance shows no general periodic behavior with slab thick- 
ness. It is of interest to note that the surface wave conductance can be appreciable for 
relatively small values of kpa, as shown in figure 3(a). For kpa = 0.595, over 90 per- 
cent of the power supplied to the aperture can be trapped within the dielectric layer. On 
the other hand, for some combinations of kpa and nonzero slab thicknesses, the surface 

wave conductance is negligible. In particular, at kpa = 1.800, there is essentially no 

z 0 

trapping in the range 0.4 < -^ < 0.6. 

A e 

The capacitance of the antenna initially rises with increasing slab thickness, oscil- 
lates, and the susceptance actually becomes inductive over a range of thicknesses for 
kpa greater than 1.305. 

The admittance for the case where the external medium is a half-space of dielec- 
tric constant 6 /€q = 2.57 (i.e., Zq = °°) is indicated in figure 3 for the first few values 
of kga. It can be seen that extremely thick slabs are required for the admittance to 
approach that of an infinite half- space. 

As a partial check, the no- coating calculations (Zq = 0) are compared with those 
of reference 1 in figure 4 for b/a = 1.57, 2.00, and 2.36. The results generally agree; 
however, there is some discrepancy which is probably due to small errors in the manual 
procedure used in reference 1 to evaluate the integrals numerically. This supposition 
seems to be supported by the more recent calculations of Curtis (ref. 5). 
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DISCUSSION 


Results have been plotted as a function of thickness in wavelengths to demonstrate 
the gross changes that can occur in the admittance after a reentry vehicle has partially 
ablated. The dielectric constant of 2.57, however, is not necessarily representative of 
that in ablation materials. This numerical value, e/6g = 2.57, was selected because it 
corresponds to that of lucite or polystyrene, both of which are commercially available in 
smooth sheets of sufficiently large area suitable for experiments at C-band and X-band 
frequencies. The loss tangents of lucite and polystyrene are also small enough at these 
frequencies to be considered negligible. 


If experimental comparisons with the computations are performed, a few comments 
are appropriate. The results of figure 3 should first be cross-plotted to give admittance 
as a function of kg a with a/Zg as a constant parameter. This should be done because 
it is much more convenient, from the experimental viewpoint, to select a constant thick- 
ness and vary the frequency. 


The dielectric constant of the material inside the line was taken as 2.00, which is 
close to that of teflon. Teflon was considered because it offers mechanical support and 
is a low-loss material which is relatively easy to machine. If other dielectrics are 
required, the results must be adjusted by multiplying each value of the admittance by 


where 


€ / e 0 


is the dielectric constant of the selected material. Care should be 


M I 
h e o’ 

exercised to keep the value of kga below the cut-off frequency of the next higher order 
mode, as determined by equation (13). 


The area of the ground plane and dielectric cover should be large enough to prevent 
the sudden termination from influencing the input admittance. The degree of influence 
of the slab and ground-plane terminations on the admittance can be ascertained by running 
a finger along the edges and noting any deflection of the VSWR meter. Circular slabs 
should be avoided because of a tendency of reflected power at the edge to be focused back 
into the aperture. 


The aperture can probably best be fed by tapering or stepping the inner and outer 
conductors of a small commercial line up to the desired size, as discussed by Ragan 
(pp. 305-314 of ref. 6). However, it is doubtful that wide-band operation can be achieved. 
It is therefore suggested that admittance measurements be taken from a slot cut into the 
large transmission line, as was done by Hartig (ref. 2). 

The variation in the characteristic impedance due to variations of the outer radius 
6b and of the inner radius 6a is given by 


6Z C 1 /5b _ 6a\ 

z c In ^ v> a / 

d 


(26) 
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If the inner radius a is 3/8 inch, then for b/a = 2.00, a variation of less than 1 per- 
cent can be expected in the characteristic impedance of the line if the tolerances in b 
and a are less than ±0.002 inch. Tolerances in concentricity can usually be neglected 
if reasonable care has been exercised in the machine work (see p. 30- 5 of ref. 7) . 

Calculations of the radiation patterns have not been included because it is doubtful 
that they could be reproduced experimentally except for those thicknesses where the sur- 
face wave conductance is negligible. If the structure is curved or terminated, that power 
which is theoretically a surface wave on an infinite structure will be released into the 
radiation field. Its interference with the computed pattern can usually be discussed only 
in a qualitative manner. 


CONCLUDING REMARKS 

The input admittance of a coaxial transmission line opening onto a flat, dielectric- 
covered ground plane is derived, and computational results are given as a function of 
slab thickness in wavelengths. A dielectric constant of 2.57 was selected for the slab, 
which physically corresponds to that for lucite or polystyrene. 

Surface wave contributions are considered in the calculations, and the results show 
that over some ranges of slab thickness 90 percent of the power supplied to the antenna 
is trapped in the dielectric. 

Comments pertinent to a physically realizable antenna are included. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., March 23, 1967, 

129-01-03-03-23. 
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APPENDIX A 


PARSEVAL'S THEOREM FOR HANKEL TRANSFORMS 


From equation (16), the complex conjugate of power radiated by the aperture can 
be given by 


P* = 7T 



E p,ap <P,°) V(P,°)P dp 


(Al) 


If equation (Al) is used directly in the computation of the aperture admittance, two 
integrations are required. First, the transform (£,0) must be integrated over all 
modes £ to obtain H^,*(p,0). Then the integration over p must be performed in 
accordance with equation (Al). Parseval's theorem reduces the problem to a single 
integration. 

To develop Parseval's theorem for Hankel transforms, first let 


E p!ap^ 0 ) = I 0 «EpJa p<€,°>Ji<«P> d * 

H cp I (p,0)=J OjJjUpJdS 


(A2) 


Then substituting equations (A2) into equation (Al) and interchanging the order of inte- 
gration gives 


P* 



d£ d^?E p I+( ^ 0) V ( ^ 0) 


pb 

\ Ji(5p)Ji(?p)p dp 

J a. 


(A3) 


— T* 

Since Ep (p, 0) is zero on the ground plane, the integration on p can be extended 
from 0 to °°. The resulting equation is (eq. (5.22) of ref. 1) 

J i(gP) J i(fo)P dP = — (A4) 

Hence, from the property of the Dirac delta function 

P* = 7T (A 5) 

Equation (A5) is Parseval's theorem for Hankel transforms. 
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APPENDIX B 


DERIVATION OF THE SURFACE WAVE CONDUCTANCE 


An inspection of the admittance integrands (eqs. (20a) to (20c)) will show that the 
denominator of equation (20b) can vanish for real fi, specifically at those values |3 = |3 n 
which are the roots of the transcendental equation 


F " V 


(Bl) 


It can further be shown that near /3 = j3 n , the denominator behaves as 0 - j3 n and 
therefore represents a pole of order one. 

For 1 S 0 = N, the integrand Fg(y3) can be expressed in the form 


F«(0) 


L(/3) /3 k € 


p 

K e cos K e + N*/c sin K e 

J 0 ( k 0 b ^) “ 

JqI 

(V^) 


[N 2 /C COS Kg 

- K e sin Ke 



(B2) 


where, for convenience, k and K e have been introduced so that 

-\ 


K = k 0 Z 0 

K e = fa ~ P\ z 0 




(B3) 


and 


p 

L(/3) = N' k cos K e - k q sin K e 

Since the poles are of order one, the integral of equation (B2) may be written 

r N V- r>N y- H/fi \ pN 

J F 2 (/3)d/3 = 7Tj 2 ResOn) + PP j F 2 (/3)d/3 = irj ^ ttj ^ + PP J F 2 (/9) d /3 

1 1 1 


(B4) 


l/3=0n 


(B5) 
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APPENDIX B 


The derivative of L(/3) at /3 = /3 n becomes 


W - 

L J ^n 


0n( k o z o) 

N 2 

cos K e,n 

% 


2 K n N 

cos K e,n + ■«— sin K e,n cos K e,n 


Sin K _ COS K 9 

+ ^ + cos 2 Ke>n 


K 


e,n 


where /c n and k q n are k and /c e , respectively, evaluated at j3 = /3 n . 
With the use of equation (Bl), equation (B6) can be reduced to 


PdL(j3)l _ M k o z o) 2 

1 ( N 2 - l\ sin 2li e,n 

L W J/3=/3 n cos x e>n 

U 2 - ij *«*» J 


In terms of K n and k q n , the numerator H^3 n ) is 

-jfk z ^ 5 

H N = ~ A° ° Kn cos K e,n + n2 % sin K e Jpo( k O a &i) - J o( k o b/ 3 n)l 
p n e,n L JL J 


And, again with the use of equation (Bl), the preceding expression reduces to 

,2 


H 


_ -3 k Q z o[ J o( k Q b M- J 0( k 0 a ^n)I 


M ~ p n cos Ke>n 

Therefore, it follows that the surface wave conductance is 

N 2 .V R fo s _ tiN 2 V [ J 0( k 0 b ^n) ~ J o( k O a/3 n)] 2 

Bs = *• - 5 " ^ es( n) ’ ' W i" 1 5 Jj + ^ n2 _ ' Sin ( 2Vo pr^| 

" 2 " ^ ( 2k o z o - Ai 2 ) _ 


(B6) 


(B7) 


(B8) 


(B9) 


(BIO) 
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